Biophysical Journal Volume 108 January 2015 251-260 251 


Spatiotemporal Detection and Analysis of Exocytosis Reveal Fusion 
“Hotspots” Organized by the Cytoskeleton in Endocrine Cells 


Tianyi Yuan,' Jingze Lu,” Jinzhong Zhang,* Yongdeng Zhang,” and Liangyi Chen'*” 

'The State Key Laboratory of Biomembrane and Membrane Biotechnology, Beijing Key Laboratory of Cardiometabolic Molecular Medicine, 
Institute of Molecular Medicine, Peking University and National Center for Nanoscience and Technology, Beijing, China; and @National 
Laboratory of Biomacromolecules, Institute of Biophysics, Chinese Academy of Sciences, Beijing, China 


ABSTRACT Total internal reflection fluorescence microscope has often been used to study the molecular mechanisms under- 
lying vesicle exocytosis. However, the spatial occurrence of the fusion events within a single cell is not frequently explored due to 
the lack of sensitive and accurate computer-assisted programs to analyze large image data sets. Here, we have developed an 
image analysis platform for the nonbiased identification of different types of vesicle fusion events with high accuracy in different 
cell types. By performing spatiotemporal analysis of stimulus-evoked exocytosis in insulin-secreting INS-1 cells, we statistically 
prove that individual vesicle fusion events are clustered at hotspots. This spatial pattern disappears upon the disruption of either 
the actin or the microtubule network; this disruption also severely inhibits evoked exocytosis. By demonstrating that newcomer 
vesicles are delivered from the cell interior to the surface membrane for exocytosis, we highlight a previously unappreciated 
mechanism in which the cytoskeleton-dependent transportation of secretory vesicles organizes exocytosis hotspots in endo- 


crine cells. 


INTRODUCTION 


In secretory cells such as neurons and endocrine cells, tran- 
sient depolarization induces Ca? entry, followed by the 
rapid fusion of secretory vesicles with the plasma mem- 
brane, thus releasing neurotransmitters and hormones to 
mediate important physiological processes (1). Electrophys- 
iological techniques, such as membrane capacitance mea- 
surements and amperometric recordings, can detect fusion 
of single vesicles with high temporal resolution (2). By us- 
ing a combination of flash photolysis, electron microscopy, 
and genetic manipulation, many aspects of the molecular 
mechanism of regulated vesicle exocytosis have been re- 
vealed (3). However, electrophysiological methods provide 
little spatial information about vesicle fusion and cannot 
observe motions of secretory vesicles before exocytosis. 
Fluorescent imaging methods can map the spatial profile 
of discrete exocytic events. Using fluorescent dyes such as 
acidic orange and FM1-43, exocytosis of acidic vesicles 
are observed in endocrine and neuronal cells (4,5). By imag- 
ing pancreatic islets in extracellular solution containing 
nonpermeable fluorescence dextrans under two-photon mi- 
croscopy, secretions buried deep within the pancreatic islets 
can be detected (6). However, the specificity of these label- 
ing protocols remains doubtful. For example, acidic orange 
has been found to localize in the acidic compartment not 
colocalized with granules (7), and extracellular labeling 
cells with fluorescence dextrans cannot distinguish between 
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exocytosis and endocytosis. Specific labeling of secretory 
vesicle exocytosis can be achieved by tagging the vesicle 
luminal cargos or vesicular membrane proteins with ge- 
netic-coded fluorescent proteins that change fluorescence 
intensity at a pH ranged from 5.5 to 7.0, such as pHluorin 
and Venus (8-10). They are quenched in the acidic vesicular 
lumen, and become dequenched and brightening in the 
neutral extracellular solution once the vesicle fusion pore 
opens, which improves the contrast of secretion signal. 
Although confocal, spinning-disc confocal, or two-photon 
microscopy can be used to detect discrete vesicle fusion 
events (11), the signal/noise ratio (SNR) of such a fluores- 
cence imaging method is compromised due to the relatively 
large excitation volume along the axial dimension. To 
further confine the focal illumination volume, total internal 
reflection fluorescence (TIRF) microscopy was developed 
(12) and used to study the dynamic behaviors of secretory 
vesicles before and during exocytosis with excellent 
contrast and better temporal resolution (4). Subsequently, 
TIRF microscopy becomes the gold standard method to 
study both regulated and constitutive vesicle exocytosis in 
a variety of cell types (13-16). 

Despite the widespread application of TIRF microscopy, 
quantitative analysis of the large amount of data generated 
by time-lapse imaging poses a challenge. It is almost impos- 
sible to manually detect and analyze the hundreds of vesicle 
fusion events recorded from single cells upon stimulation un- 
der a TIRF microscope. Most researchers rely on the manual 
annotation of a limited number of fusion events. Such anal- 
ysis is prone to the biases of selection and does not always 
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lead to a statistically supported conclusion. Recently, a few 
groups have started to develop algorithms that facilitate the 
identification of vesicle fusion from time-lapse images. For 
example, Bai et al. and Huang et al. reported programs that 
enable direct analysis of the docking and fusion kinetics of 
glucose transporter 4 (GLUT4) storage vesicles (GSVs) 
(13,17). However, these methods are semiautomatic and 
require extensively manual inspection and revision of indi- 
vidual events. Sebastian et al. (18) implemented an auto- 
mated algorithm that extracts the spatial location and onset 
time of each fusion by a forward subtraction method. Such 
an algorithm does not fully use the time-sequential informa- 
tion from image stacks. Therefore, although it could detect 
~86% of the true fusion events, the specificity was only 
65%. Based on particles tracking and statistical testing of 
the similarity between candidate events and true fusion 
events, two other algorithms were proposed, but the rate of 
false positive events was even higher with noisy images 
(19,20). Hence, none of these methods is widely used. 
Furthermore, except for one (18), none of these works take 
full advantage of the spatial information available to conduct 
spatial analysis of all vesicle fusion events. 

The release of synaptic vesicles in synaptic transmission is 
spatially confined to presynaptic terminals. Abundant synap- 
tic vesicles cluster at the densely packed presynaptic region 
(active zone), which is organized around scaffolding proteins, 
such as ELKS and Rab3-interacting molecule (RIM), and 
these proteins contribute to the spatial preference (21). Iso- 
forms of these proteins also exist in endocrine cells such as 
pancreatic 6-cells. Although no active zone-like physical 
structures have been observed in these cells under electron mi- 
croscopy, secretory vesicles have been found to sequentially 
fuse at the same spot, and this behavior is proposed to depend 
on ELKS proteins (22). Nevertheless, whether vesicles 
repeatedly fuse at preferred locations on the plasma mem- 
brane has not been statistically proven. Cytoskeletons may 
play a role in the spatial coordination of the release of secre- 
tory vesicles. By acting either as tracks for the transportation 
of secretory vesicles (23) or as barriers that prevent collapse of 
vesicles with the plasma membrane (24), cortical actin has 
been proposed to facilitate or inhibit regulated exocytosis. 
Similarly, microtubule cytoskeleton may also participate in 
secretion processes. Cortical microtubules extend within 
100 nm underneath the plasma membrane, and vesicles 
move along microtubules beneath the plasma membrane until 
fusing with the plasma membrane (25). Consistently, sponta- 
neous neurotransmitter release in excitatory synapse was 
increased or inhibited in neuron overexpressing or knock- 
down of cytoplasmic-linker-associated protein2 (CLASP2), 
a microtubule-associated protein that targets microtubules 
to the cell cortex (26). Despite these experiments, how disrup- 
tion of cytoskeleton network quantitatively affects the spatial 
preference of fusion events remained to be explored. 

Here, by combining two different fusion detection para- 
digms, we have developed an algorithm to detect vesicle 
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fusion based on a sudden increase in the fluorescence signal 
after a fusion pore opens. We have demonstrated the supe- 
rior sensitivity and accuracy of our algorithm in detecting 
vesicle fusion from different cell types. Using this platform, 
we revealed a tightly orchestrated cortical cytoskeleton 
network that delivers secretory vesicles for fusion at 
spatially confined regions in endocrine cells. 


METHODS 
Materials 


VAMP2-pHluorin, IRAP-pHluorin, and VAMP2-EGFP plasmids were 
kindly provided by Prof. Tao Xu. INS-1 and 3T3-L1 cells were cultured 
as previously described (27,28). Cells were transfected using Lipofectamine 
2000 reagent (Invitrogen, Carlsbad, CA) according to the manufacturer’s 
instructions. Cells were detached using trypsin-EDTA 24 h posttransfection, 
plated onto poly-L-lysine-coated coverslips, and cultured for 20-28 h before 
experiments in a 37°C, 5% CO, incubator. Cytochalasin D (C8273) and 
nocodazole (M1404) were from Sigma-Aldrich (St. Louis, MO). 


Live-cell imaging 


All experiments were performed at 37°C. Before an imaging experiment, 
INS-1 cells were washed and incubated on glass coverslips in a bath solu- 
tion containing (in mM) 136 NaCl, 4.2 KCl, 2.4 CaCl, 1.2 KH2PO,, 1.2 
MgSO,, 4 glucose, 10 HEPES, and 1 L-glutamine (pH 7.4), and the cover- 
slip was placed in a metal chamber mounted on a heated stage. Throughout 
the study, INS-1 cells were stimulated with a high glucose and high K* so- 
lution containing (in mM) 70 NaCl, 70 KCl, 2.4 CaCl», 1.2 KH2PO,, 1.2 
MgSO,, 15 glucose, 10 HEPES, and 1 L-glutamine (pH 7.4) to trigger 
exocytosis. 3T3-L1 cells were kept in KRBB solution containing (in 
mM) 129 NaCl, 4.7 KCl, 1.2 KH,PO,4, 5 NaHCOs, 10 HEPES, 3 glucose, 
2.5 CaCl», 1.2 MgCl, and 0.1% bovine serum albumin (pH 7.2); 100 nM 
insulin was applied later to trigger exocytosis (13). 

For actin filament disruption, 20 uM cytochalasin D (CytoD) was added 
to the culture medium for 20 min before the imaging experiments. For 
microtubule disruption, 10 uM nocodazole (Nocod) was added to the cul- 
ture medium for 1 h before the imaging experiments. 


The TIRF microscope 


For Figs. 1—4, a TIRF microscope equipped with a Sensicam EM682KX 
(PCO AG, Kelheim, Germany) and a 100x, 1.49 numerical aperture, oil- 
immersion TIRF objective (Olympus) was used. Images were acquired at 
2 Hz. Each pixel was 67 x 67 nm. For Fig. 5, a different TIRF setup was 
used. We used an Olympus IX81 inverted microscope equipped with a 
150 x, 1.45 numerical aperture, oil-immersion TIRF objective (Olympus). 
The fluorescence signals were acquired by an electron-multiplying charge- 
coupled device camera (Andor iXon3 888) at a sampling rate of 6 Hz 
that was controlled using MetaMorph (Molecular Devices). Each pixel 
from the Andor camera was 87 x 87 nm. A 473 nm laser was used to excite 
EGFP/pHluorin. 


RESULTS 


Automatic identification and analysis of fusion 
events 


Fusion of fluorescently labeled vesicles with the plasma 
membrane observed with TIRF microscopy has been 
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described as a sudden brightening and rapid spreading of 
the fluorescence signal (29). The sudden increase in the 
fluorescence signal occurs because evanescent excitation 
is more efficient closer to the interface and because of 
the dequenching of the pHluorin fluorescence once the 
acidic granule lumen has been exposed to a neutral envi- 
ronment when the fusion pore opens (9). The sudden in- 
crease in the fluorescence signal is followed by an 
exponential decrease in fluorescence intensity, which is 
the result of the rapid diffusion of vesicular membrane 
proteins or contents out of exocytic regions. We used a 
combination of these features (fluorescence increase, 
decrease, and ensuing diffusion of fluorophore) to detect 
vesicle fusion. The framework of the algorithm consisted 
of five major steps (Fig. | A). First, the image sequence 
was denoised, and the background was subtracted using 
the a trous wavelet transform algorithm (30) (step 1). 
Here, a hard threshold parameter (the selected wavelet 
component) after wavelet transformation ultimately deter- 
mines the sensitivity and the accuracy of the spot detection 
algorithm. Therefore, when the imaging condition changes, 
we always manually adjust this parameter at the first 
image from a series of image stacks to ensure that proper 
numbers of fluorescence spots were detected. Once deter- 
mined, the parameter can be used for the rest of the images 
with similar quality without much adjustment. The locali- 
zation of all fluorescent particles (representative of vesi- 
cles) was then extracted with either a local maxima or a 
weighted centroid method (within a range of 5 x 5 pixels) 
at each frame (step 2). After particle detection, time-lapse 
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trajectories of these vesicles were constructed using a mul- 
tiple particle tracking software based on the linear assign- 
ment problem algorithm (31) (step 3). Here, our software 
pops out four parameters that may need to be adjusted ac- 
cording to vesicle behavior, including the gap closing fac- 
tor, the minimum and maximum search radius, and the 
particle motion type (linear or not). Zero gap closing will 
increase the true positive rate (TPR) of the fusion detector; 
conversely, increasing the gap closing factor reduces the 
false positive rate (FPR) at the expense of sensitivity. 
The start and the end of a trajectory represents the appear- 
ance and the disappearance of a vesicle. Because vesicle 
fusion led to the end of a trajectory, we probed a possible 
vesicle fusion event within one track (step 4), which 
considerably reduced the time required for fusion detection 
(19,20). This step also helps to eliminate false positive 
fusion detection due to the lateral movement of vesicles. 
Finally, we performed a two-dimensional Gaussian fit 
and spatiotemporal analysis to extract additional informa- 
tion (step 5). 

Fusion detection is the most difficult step and we used a 
protocol that combined advantages of previous algorithms. 
First, we selected vesicles that once exhibited fluorescence 
intensity (F,) that was greater than the background fluores- 
cence (Fo) by a certain ratio (R) within its lifetime as 
possible fusion events (Fig. | B). The background was 
calculated as the mean intensity of the area (within a range 
of 5 x 5 pixels) at the center of the punctum (probable 
vesicle) and was measured for >5 frames before the appear- 
ance or >5 frames after the disappearance of the punctum. 


mm VAMP2 

IRAP 

FIGURE 1 Automatic identification and anal- 
ysis of fusion events. (A) The framework of the 
algorithm, comprising five steps (non-rounded 
rectangles). (B) An example of a typical fusion in 
an INS-1 cell labeled by VAMP2-pHluorin (lef 
and its central intensity profile (right). The rapid 
increase in the central intensity indicates the open- 
ing of a fusion pore; thus, the peak/background ra- 
tio (parameter R,) can be used to detect candidate 
events. (C) Consecutive image subtraction of (B) 
and its central intensity profile; only those peak in- 
creases that are greater than the background fluctu- 
ation by R, will be considered candidate events. 
The montages are 19 x 19 pixels shown at 0.5-s in- 
tervals, and the corresponding kymographs (below) 
are 150 s. (D) The performance of the method for 
pHluorin-labeled fusion events (n = 1701 events 
from 6 cells for VAMP2 and n = 531 events 
from 7 cells for IRAP). R; and R, were set to 1.3 
and 3 in both situations. Error bars are + SE. 
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The method was basically a time-lapse one-dimensional in- 
tensity analysis (19,20), prone to error when local SNR near 
fusion sites decreased after extensive exocytosis. In addi- 
tion, we used consecutive image subtraction to normalize 
the background fluorescence intensity across the entire field 
of view (18), and detected possible fusion events as a local 
increase in the fluorescence intensity over a preset ratio 
(R2) (Fig. 1 C). The detection accuracy and efficiency 
were much better with a combination of both R; and R3 
than using only one of these parameters (Fig. S1 in the 
Supporting Material). 

In insulin-secreting cells, insulin is stored in secretory 
granules and is released upon an elevation in glucose levels. 
We used INS-1 cells containing secretory vesicles labeled 
with VAMP2-pHluorin, and detected massive vesicle fusion 
triggered by stimulation with 15 mM glucose and 70 mM 
KCl. By tuning R; and R3 to proper values (1.3 and 3 in 
this set of data), we detected vesicle fusion events with a 
TPR of 96.5% and a FPR of 9.5% (Fig. 1 D); this FPR is 
much lower than that of previous methods (18-20). Next, 
we tested whether our algorithm could be used to detect 
other types of vesicle fusion. Insulin released into the blood- 
stream acts on the insulin receptors on adipocytes and mus- 
cle cells. This triggers the GSVs to fuse with the plasma 
membrane, leading to the insertion of the GLUT4 protein 
into the plasma membrane and potentiating the glucose 
uptake by these cells (32). Unlike Ca**-dependent insulin 
secretion, GSV fusion depends on an insulin-stimulated 
signaling pathway but not on the cytosolic Ca?™ concentra- 
tion. We labeled GSVs with pHluorin that was tagged to the 
luminal side of insulin-regulated aminopeptidase (IRAP, 
IRAP-pHluorin), a transmembrane protein that colocalizes 
with GLUT4 on GSVs. In adipocyte-like 3T3-L1 cells, insu- 
lin stimulation (100 nM) slowly triggers the fusion of GSVs 
(33). We still achieved a TPR of 93.5% for GSV fusion 
(Fig. | D). However, the FPR from the IRAP-pHluorin ex- 
periments was higher than that from the VAMP2-pHluorin 
experiments, possibly due to the inferior image quality 
with high levels of noise and background fluorescence. 
When setting the parameters, there was always a compro- 
mise between TPR and FPR. With high FPR, we could in- 
crease the values of R; and R, to eliminate false positive 
nonfusion events. With low TPR, we selected a fusion event 
that was missed in previous detection because of low local 
SNR, and the program would automatically update current 
R; and R, to the lowest values between values of this fusion 
event and preexisting values. These updated parameters 
could be used to detect fusion events from the same image 
stack thereafter. By optimizing in this manner for a few 
rounds, it is easy for the researcher to reach the best compro- 
mise between TPR and FPR for the image stack tested. 
The other parameters in the software package (including 
denoising and tracking) would not vary much for images 
acquired under similar conditions. Subsequently, the rest 
of the images from the same batch of experiments could 
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be processed with little adjustment in a batch-processed 
manner with high sensitivity and accuracy. This procedure 
ensures our program to be adaptive and robust in detecting 
different types of exocytic events from different cells. 


Disruption of the cytoskeleton inhibits regulated 
exocytosis and changes the fusion kinetics 


Next, we evaluated whether our software could be used to 
quantitatively study exocytosis in INS-1 cells. When vesi- 
cles were labeled with VAMP2-pHluorin, a typical fusion 
event exhibited an abrupt increase in fluorescence followed 
by a rapid decrease, which was fitted with an exponential 
function (Fig. 2 A). In six cells stimulated with high concen- 
trations of KCl and glucose, we identified up to ~1700 
fusion events. Using the results, we assembled the average 
time course of the change in exocytosis after stimulation 
(Fig. 2 B), the average fusion events per um? per s (Fig. 2 
C) and the average decay in fluorescence after fusion 
(Fig. 2 D). Disruption of either the actin filaments using Cy- 
toD or the microtubule network using Nocod significantly 
inhibited the time-dependent fusion of vesicles upon stimu- 
lation (Fig. 2 B). Overall, average fusion rates were reduced 
to one-third in the CytoD-treated cells and one-sixth in the 
Nocod-treated cells, respectively (Fig. 2 C), indicating that 
both actin and microtubules play an indispensable role in 
transporting secretory vesicles to the cell surface before 
exocytosis. Unexpectedly, both the CytoD and Nocod treat- 
ments also caused a dramatic increase in the fusion decay 
time (Fig. 2 D), suggesting that the cytoskeleton affects 
either the fusion pore dilation or the protein diffusion on 
the plasma membrane. After performing a two-dimensional 
Gaussian fit of the fluorescence intensities at different times 
after the fusion pore opened (Fig. 2 E), we directly plotted 
the squared widths of these Gaussian functions versus 
time and determined the apparent diffusion coefficient of 
VAMP2-pHluorin on the plasma membrane (34). Consistent 
with the change in the decay time, the apparent diffusion co- 
efficient decreased in both the CytoD- and Nocod-treated 
cells (Figs. 2 F). Therefore, disruption of the cytoskeleton 
networks slowed down the diffusion of the vesicle mem- 
brane proteins on the plasma membrane after the fusion 
pore opened. 


The cytoskeleton is essential for spatially 
clustered exocytosis in INS-1 cells 


In addition to the temporal occurrence, a detailed spatial 
map of exocytosis initiation is also provided from these im- 
ages. It has been reported that some vesicle fusion sites were 
more likely to harbor another fusion event within a granule 
diameter (~300 nm) in chromaffin cells (29). The phenom- 
enon, called sequential exocytosis, has been found in many 
other cells (35-39). Therefore, we studied the spatial distri- 
bution of fusion events in INS-1 cells. It was apparent that 
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some vesicles fused with the plasma membrane at similar 
locations (within a 5 x 5 pixel area) at different times 
(sequential exocytosis, Fig. 3 A). The sequential fusions 
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FIGURE 2 Disruption of the cytoskeleton in- 
hibits secretion and changes the fusion kinetics. 
(A) An example of a fusion event and the time- 
dependent fluorescence intensity profile. The mon- 
tages are 19 x 19 pixels at time intervals of 0.5 s, 
and the corresponding kymograph (below) is 150 s. 
The decay time 7 of the central intensity was ob- 
tained by fitting with an exponential function. (B) 
Time-dependent changes in the exocytosis rate in 
control cells and in cells pretreated with either Cy- 
toD or Nocod. A KCI and glucose solution was 
applied to the tested cells at time 0. (C) Average 
fusion rate for each condition. (n = 1701 events 
from 6 cells for control, n = 824 events from 6 cells 
for CytoD treatment, and n = 577 events from 7 
cells for Nocod treatment). (D) Average decay 
time for each condition. (n = 911 events from 6 
cells for control, n = 622 events from 6 cells for 
CytoD treatment, and n = 331 events from 7 cells 
for Nocod treatment). (E) Intensity profiles of an 
example fusion event at different times fitted by 
Gaussian functions. (F) Average apparent diffusion 
coefficients for each condition. (n = 1629 events 
from 6 cells for control, n = 826 events from 6 cells 
for CytoD treatment, and n = 467 events from 7 
cells for Nocod treatment). Error bars are + SE, 
and all statistical significance values are relative 
to the control. p values were determined using 
Student’s t-test. ** p < 0.01; *** p < 0.001. 


constituted 28% of all fusion events in control cells, but 
this proportion dramatically decreased to 10% of the events 
in cells treated with either CytoD or Nocod (Fig. 3 B). In 


FIGURE 3 Disruption of the actin or microtu- 
bule network reduces the ratio of sequential fu- 
sions. (A) An example of sequential fusion and 
the central intensity profile. The interval between 
fusions was calculated as AT. The montages are 
19 x 19 pixels at time intervals of 0.5 s, and the 
corresponding kymograph (below) is 150 s. (B) 
The ratio of sequential fusion for each condition. 
(C) Average minimal distance between each fusion 
event. In control cells, the minimum distance be- 
tween each event was 10 pixels (n = 1701 events 
from 6 cells for control, n = 824 events from 6 cells 
for CytoD treatment, and n = 577 events from 7 
cells for Nocod treatment). (D) Average AT for 
each condition (in B and D, n = 478 events from 
6 cells for control, n = 101 events from 6 cells 
for CytoD treatment, n = 65 events from 7 cells 
for Nocod treatment). Error bars are +SE. 
p values were determined using Student’s t-test. 
*** p < 0.001. 
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agreement with the selective reduction of repeated exocy- 
tosis at the same location, the average distance between 
two nearest fusion events increased as the actin or microtu- 
bules were disrupted (Fig. 3 C). The average time interval 
between two sequential fusions was longer in either the Cy- 
toD- or Nocod-treated cells compared with that in the con- 
trol (Fig. 3 D), reflecting a reduced fusion efficiency even 
for vesicles undergoing sequential exocytosis. Overall, 
vesicle fusion becomes more sparsely distributed when the 
cellular cytoskeleton is disrupted. 

Even for the vesicles undergoing sequential exocytosis, 
the spatial localizations of different fusion sites did not al- 
ways precisely overlap. In this regard, we wonder whether 
randomly diffusing vesicles could hit and fuse at the same 
area in the plasma membrane multiple times by chance and 
contribute to the sequential fusion. To address this ques- 
tion, we used Ripley’s K-function (40) to statistically 
examine the spatial distribution of all fusion events within 
a cell (Fig. 4 A). The spatial K-function was determined by 
counting the number of fusion events within a specific dis- 
tance from one event (18). From 100 Monte Carlo simula- 
tions using the same number of fusion events randomly 
distributed within the cell boundary, we constructed a 
range of K-function variations at specific distances. The 
K-function trace of the fusion locations from the real 
experiment were well beyond that range, indicating that 
the exocytosis was not randomly distributed in INS-1 cells 
(Fig. 4 B). Disruption of either actin or microtubules 
eliminated the clustering of fusion sites and spatially ran- 
domized the fusion events (Fig. S2, A and B). Next, we 
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introduced a parameter, the uniform index, to quantitatively 
evaluate the degree of spatial aggregation of exocytosis. By 
creating a Voronoi diagram to separate each fusion event, 
we determined the average area occupied by each fusion 
(Fig. 4 C and Fig. S2 C). The CytoD or Nocod treatment 
increased the average area occupied by each fusion, in 
agreement with the severe inhibition of sequential exocy- 
tosis. The uniform index was calculated as the standard de- 
viation of the fusion area divided by the mean of the fusion 
area. Monte Carlo simulation experiments were performed 
by randomly spreading the fusion events over the cell sur- 
face area using the number of fusion events within the 
same cell boundary (Fig. S2 D). Compared with the simu- 
lation, the vesicle fusion in real experiments exhibited a 
63% higher uniform index in the control cells but was 
not significantly different in the cells treated with CytoD 
or Nocod (Fig. 4 D). Considering these findings, we 
conclude that the vesicle exocytosis occurs at spatially 
preferred regions in INS-1 cells (hotspots) and that this 
preference is dependent on the presence of an intact 
cellular cytoskeleton network. 


Disruption of the cytoskeleton selectively 
reduces the fusion of newcomer vesicles 


Despite the establishment of the presence of cytoskeleton- 
dependent exocytosis hotspots in INS-1 cells, it was unclear 
how the cytoskeleton mediated this process. To address this 
question, we changed the vesicle marker from VAMP2- 
pHluorin to VAMP2-EGFP, which is used to observe the 


z FIGURE 4 Actin and microtubule networks 
ree regulate the spatial pattern of cell secretion. (A) 
A control cell with all fusion events from a time- 
lapse series marked on the first frame (crosses). 
Scale bar, 5 um. (B) Ripley’s K function of (A). 
Solid line, observed spatial K-function of fusion 
events; dashed lines, upper and lower envelopes 
of 100 simulation results with the same number 
of random events within the cell boundary as in 
(A). (C) Average area occupied by each fusion for 
each condition. (D) Average uniform index for 
each condition (in C and D, n = 1701 events 
from 6 cells for control, n = 824 events from 6 cells 
for CytoD treatment, and n = 577 events from 7 
cells for Nocod treatment). Error bars are +SE. p 
values were determined using Student’s t-test. 
RK  < 0.001. 
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trafficking, tethering, docking, and priming of vesicles 
before fusion. Using TIRF microscopy, we found that fusion 
events can be classified into three distinct categories (41) 
(Fig. 5, A-D). Some vesicles were already docked on the 
plasma membrane for a period of time before stimulation 
and then fused, called predock fusions (Fig. 5 A). Some ves- 
icles newly entered the evanescent field and fused with the 
plasma membrane immediately without evident docking, 
called newcomer-no docking (Fig. 5 B). Some vesicles 
newly entered the TIRF zone and underwent a short period 
of docking and then fused, called newcomer-short docking 
(Fig. 5 C). Disruption of the cytoskeleton did not inhibit 
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the predock fusions but significantly inhibited the fusion 
rates of newcomer-no docking and newcomer-short docking 
(Fig. 5 E). In addition, the density of stably docked vesicles 
did not change after cytoskeleton disruption (Fig. S3). 
Therefore, rather than mediating the accumulation of secre- 
tory vesicles under the plasma membrane, the cytoskeleton 
network delivers secretory vesicles from the cell interior to 
specific regions on the cell surface membrane. 

Next, we analyzed the docking latency in newcomer-short 
docking events by fitting the intensity profile before fusion 
with a step-finding approach (Fig. 5 D). The docking time 
was shortened in the CytoD-treated cells, consistent with 
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dock fusion. (D) Central intensity profiles of (A— 
C). For newcomer-short dock fusion, the interval 
between the two dashed lines in (D) was calculated 
as the docking time. (E) Average fusion rate of 
each type of fusion for each condition. (n = 1016 
events from 6 cells for control, n = 476 events 
from 9 cells for CytoD treatment, and n = 285 
events from 8 cells for Nocod treatment). (F) 
Average docking time for each condition (n 
131 events from 3 cells for control, n = 150 events 
from 9 cells for CytoD treatment, and n = 79 
events from 8 cells for Nocod treatment). (G) 
Example of lateral shift of a fusing vesicle. 
The displacement from the last position before 
fusion (black circle) to the final fusion site (white 
circle) was 114 nm. Scale bar, 500 nm. (H) Distri- 
bution of the lateral shift in control cells. (D) 
Average lateral shift for each condition. (n = 262 
events from 3 cells for control, n = 223 events 
from 9 cells for CytoD treatment, and n = 129 
events from 8 cells for Nocod treatment). 
Error bars are + SE. p values were determined us- 
ing Student’s t-test. * p < 0.05, ** p < 0.01, and 
*** p < 0.001. 
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previous works (25), but not in the Nocod-treated cells 
(Fig. 5 F). This finding suggests that cortical actin filaments 
prolong the physical docking process of secretory vesicles. 

Vesicles can undergo a lateral displacement under the 
plasma membrane before fusion in chromaffin cells 
(42,43). We also found such behavior in INS-1 cells, as 
shown in an example in Fig. 5 G, in which the distance 
from the position before fusion to the final fusion site was 
114 nm. Because our method contains a robust multiple par- 
ticle tracking algorithm, we were able to detect and analyze 
such behaviors systematically. The lateral shift varied from 
tens to hundreds of nanometers with a mean value of 131 + 
5.8 nm in control cells (Fig. 5 H). However, the value was 
reduced to 99 + 5.2 nm in CytoD-treated and 91 + 
6.4 nm in Nocod-treated cells (Fig. 5 1). This finding again 
supported a role of the cortical cytoskeleton in arranging the 
final spatial redistribution of vesicle fusion events. 


DISCUSSION 


Here, we have developed a computer-assisted identification 
and analysis software to detect pHluorin- and EGFP-labeled 
fusion events in INS-1 cells and 3T3-L1 cells, which pro- 
vides superior sensitivity and accuracy as compared to pre- 
vious methods (18—20). Our software can also be used to 
detect vesicle fusion events labeled by other markers, such 
as VGLUT1-mOrange2 (44), GFP-insulin (45), and NPY- 
pHluorin (46), which all share similar fluorescence charac- 
teristics during fusion. Prolonged acquisition always poses a 
challenge to automated image processing, as morphologic 
changes in the cell shape might result in changes in focus 
and reduction in image contrast. Our algorithm is capable 
of analyzing image stacks with a long acquisition time 
with relative robustness (Fig. S4). Because the algorithm 
is based on trajectories construction, it will effectively 
reduce false positive events exhibiting as transient, local 
disappearance of fluorescence intensity, which are possibly 
a result of the changes of focus or fluorophore blinking. The 
increase in FPR during prolonged experiments was mainly 
due to the accumulation of VAMP2-pHluorin on the plasma 
membrane after extensive exocytosis, which increased the 
cellular background fluorescence and reduced the image 
contrast during exocytosis (Fig. S4). This could be resolved 
by targeting pHluorin to luminal cargos that diffuse into the 
extracellular solution after exocytosis. Overall, after simple 
adjustment, our program is robust and adaptive to detect 
different types of vesicle fusion under different experi- 
mental conditions. In addition, the tracking algorithm can 
be used to track endocytic markers, such as clathrin and 
dynamin from time-lapse images (31,47). Therefore, it is 
also possible to modify our algorithm to study the coupling 
mechanisms between exocytosis and endocytosis in the 
future. 

Using this tool, we statistically prove that stimulated 
INS-1 cells release secretory vesicles at spatially preferred 
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hotspots. In principal, our data agree with sequential 
exocytosis found in chromaffin cells (29), pancreatic acinar 
cells (48), and islets (39). Different ratios of sequential 
exocytosis to the total exocytosis may be due to different 
cell types and preparation methods. Different microscopic 
methods may also contribute to the differences. For 
example, two-photon microscopy data were collected at a 
sampling frequency of 0.3~2 s per frame (6,39), which is 
slower as compared to the sampling speed of current study 
and may miss some of the fast fusion events. Distinct from 
the active zone in synapses, the exocytosis hotspots in 
these endocrine cells are not formed by the fusion of 
densely packed vesicles that are predocked under the 
plasma membrane. In contrast, secretory vesicles that 
fuse at hotspots are delivered to the plasma membrane 
from the cell interior via microtubular and actin filaments. 
Disruption of either the actin or microtubular network 
severely reduced the overall exocytosis in addition to 
randomizing the vesicle fusion sites, suggesting that 
spatially confined exocytosis plays a major role in stimula- 
tion-triggered hormone release from insulin-secreting cells. 
The molecular mechanism underlying the cytoskeleton- 
mediated targeted delivery of vesicles remains unknown. 
Microtubules interact with the active zone protein ELKS 
through its binding protein CLASP (49) and to recruit 
vesicles to specific membrane locations, such as focal ad- 
hesions (50). Interactions between actin and vesicle mem- 
brane proteins, such as synapsin, have also been proposed 
(51). Future investigations will determine whether one or 
all of these mechanisms participate in the hotspot secretion 
of endocrine cells. 

It was interesting to find that the disruption of the cyto- 
skeleton changed the lateral shift of vesicles during the 
final fusion process (Fig. 5 J). This result showed the exis- 
tence of a cortical cytoskeleton network and its role in 
precisely interacting with components on the plasma mem- 
brane. In agreement with this hypothesis, the diffusion of 
VAMP2 on the plasma membrane was slowed down in 
the absence of cytoskeleton filaments (Fig. 2, D and F). 
This was unlikely to be caused by the facilitatory effects 
of the cytoskeleton on fusion pore dilation because actin 
has been proposed to negatively regulate fusion pore 
expansion (52) and stabilize Q-shaped structures during 
vesicle fusion (48). Alternatively, cytoskeleton filaments 
may interact with lipid domains on the plasma membrane 
to facilitate the diffusion of membrane proteins on the 
plasma membrane (47). 

In summary, we have developed a software platform for 
the identification and spatiotemporal analysis of vesicle 
fusion events obtained by TIRF microscopy. Using this 
tool, we revealed that newcomer vesicles are delivered 
from the cell interior to the plasma membrane via cytoskel- 
eton filaments to fuse at hotspots, and these vesicles make a 
major contribution to insulin-secretion from INS-1 cells 
stimulated with glucose and depolarization. 
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